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Abstract 

An algorithm is presented here to get more detailed information, of mixed cul- 
ture type, based exclusively on the biomass concentration data for fermentation 
processes. The analysis is performed with only the on-line measurements of the re- 
dox potential being available. It is a two-step procedure which includes an Artificial 
Neural Network (ANN) that relates the redox potential to the biomass concen- 
trations in the first step. Next, a multifractal wavelet analysis is performed using 
the biomass estimates of the process. In this context, our results show that the 
redox potential is a valuable indicator of microorganism metabolic activity during 
the spontaneous fermentation. In this paper, the detailed design of the multifractal 
wavelet analysis is presented, as well as its direct experimental application at the 
laboratory level. 
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1 Introduction 



During spontaneous fermentation processes, microorganisms employ sugars 
and other constituents of must as substrate for their growth, converting them 
into ethanol, carbon dioxide, higher alcohols and their esters, and other metabolic 
compounds. In general, it could be useful to know more about the dynamics 
of the entire microflora during spontaneous fermentations. In particular, an 
analysis that enables a monitoring process could be fundamental for a quality 
control that ensures at least the homogeneity of the final product. 

Nevertheless, a bottleneck in all biochemical monitoring process is often the 
lack of sensors for biological variables. Moreover, it is a well-known issue that 
in order to monitor many biotechnological processes, the problem of growth 
rate estimations represent a strategic feature. Previously, various attempts of 
relating the redox potential to fermentation processes have been made taking 
into account that redox potential assesses the growth ability of microorganisms, 
as well as the physiological activity in a given environment [2], [3], [I], [5]. 
Particularly, the practical significance of redox potential and oxygen content at 
various stages of winemaking was examined in [BJ. Many chemical, enzymatic 
and biological processes in wine are correlated with the oxidative state of the 
wine. 

It is well known that many natural systems exhibit complex dynamics de- 
scribed by long-range power laws. In some cases the output of such systems, 
which corresponds to a fluctuating time series, may be characterized or quan- 
tified by a spectrum of exponents called the multifractal spectrum. The mul- 
tifractal formalism to characterize processes, measures, and functions, intro- 
duced by Halsey and collaborators [7j, opened a new direction for the search 
of dimension type characteristics of dynamical systems in the form of spec- 
tra for dimensions that reflect both the geometric and dynamical structure 
of nonlinear systems. Despite the fact that there are a variety of methods 
to calculate the singularity spectrum of a multifractal structure, the method 
based upon the wavelet theory appears to be the most appealing and success- 
ful [8]. In fact, the interest in wavelet methods is that they are numerically 
more stable [9]. In 1993, Bacry et al. [TU] developed this method based on the 
definition of a partition function in which the concept of wavelet transform 
modulus maxima (WTMM) was implemented. They demonstrated that the 
singularity spectrum for a Bernoulli measure or a fractal distribution can be 
readily determined from the scaling behavior of such a partition function and 
similar results have been proved for more general measures by Murguia and 
Unas p]. 

The paper is organized as follows. Section[2]is devoted to a concise presentation 
of the fermentation experiment performed to illustrate our approach. The data 
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analysis, which is a combination of an ANN and a multifractal wavelet scheme, 
is described in Section [3j Finally, the paper ends up with some concluding 
remarks. 



2 Experimental setup 

2.1 Microorganism and culture conditions 

In order to evaluate experimentally how suitable the estimation algorithm 
is, we performed six individual batch experiments with the must of Agave (or 
Agave syrup) for which we used inocula of native microorganisms (without the 
addition of any commercial strain). This must was centrifuged at 8000 rpm x 
10 min and stored in a frozen state at — 20°C prior to experiments. 

The batch fermentations were carried out in a bioreactor with pH and redox 
sterilizable electrodes. The electrodes are connected to a console for data ac- 
quisition, a device which is connected to a computer where the data are stored 
and computed as can be seen in Fig. [TJ 

The bioreactor was filled with 900 mi of must as a culture medium, 100 mi 
of the inoculum in its exponential growing phase (biomass 0.1 g/L) and 0.1 
% of ammonium sulfate at final concentration. The initial conditions of the 
fermentation were settled at a temperature of 32.5°C and initial sugar con- 
centration of 70 g/L. The pH does not show a dynamic evolution, maintaining 
itself at a value of 4 during the whole process. 



2.2 Analytical procedures and measurements 

The batch processes have been monitored for 14 hours, through sampling un- 
der sterile conditions. In order to quantify biomass and ethanol concentrations, 
5 mL samples of culture were removed every 30 min. The samples were cleared 
by centrifugation at 6000 rpm for 5 minutes at room temperature. The next 
step was to collect the supernatant phase and store it frozen at — 20°C prior to 
be analyzed. The obtained pellet was resuspended in distilled water in order 
to proceed with the biomass analysis. 

The biomass measurements have been performed using (varian) UV spec- 
troscopy at 600 nm. The obtained values were interpolated with a standard 
curve of cell dry weight concentration. 
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Fig. 1. The schematic representation of the experiments carried out in our labora- 
tory. 



2.3 Redox potential 



The measurement of redox potential is relatively fast, accurate, and reliable 
and its values give an insight into the oxidation process as well as the inherent 
ability of reduction in the process, which is well established in the case of wine 
[BJ. The measured values of the redox potential can give information on redox 
reactions in wine, which have an important effect on its quality and stability 
[BJ. During storing and aging of wine, oxidation and reduction processes affect 
the character and taste of wine to a considerable extent [BJ, [5]. In our Agave 
case, the redox potential measurements were acquired periodically each 0.01 hr 
during 14 hrs, and the data were stored and processed on line by means of a 
computer (see Fig. [Q. 



3 Data analysis 



3.1 Artificial neural network 



In order to relate the redox measurements to the biomass concentration an 
ANN procedure is applied. The methodology that we carried out includes a 
forward-propagation training algorithm for the ANN using some of our ex- 
perimental data. In order to perform our task we construct a model of the 
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following form: 



X 1 = g(X 2 ) (1) 

where X2 represents the redox potential measurement data (mV), X\ = F(t) is 
the set of time-dependent biomass concentration data (mg/L) and the function 
g(X2) is approximated by means of the ANN procedure. The ANN architecture 
is of the standard type [12] with a single ANN hidden-layer containing 10 units. 



Each unit of this network uses a sigmoid function as the activation function. 
On the other hand, the output contains a linear activation function, in our 
case the identity. The feed-forward training algorithm considered here is the 
conjugate gradient method [13j. Three of the six individual batch experiments 
were used to provide data for the training process. The ANN after the training 
gives an error of only 0.0029. 



3.2 Multifractal wavelet analysis 



The wavelet transform (WT) of a distribution function s(t) is given by 

W s (a,b) = - ( 8 (t)$( — )dt, (2) 



a J \ a 

-00 



where ip is the analyzing wavelet, & 6 R is a translation parameter, whereas 
a E R + (a 7^ 0) is a dilation or scale parameter, and the bar symbol denotes 
complex conjugation. 

One fundamental property that we require in order to analyze singular behav- 
ior is that if)(t) has enough vanishing moments [TUfTl] . A wavelet is said to 
have n vanishing moments if and only if it satisfies 

00 

J t k i()(t)dt = 0, for k = 0,1,..., n- 1, (3) 

—00 

and 

00 

J t k ip(t)dt ^ 0, for k = n. (4) 

—00 

This means that a wavelet with n vanishing moments is orthogonal to all 
polynomials up to order n—1. So the wavelet transform of s(t) with a wavelet 
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ip(t) with n vanishing moments is nothing but a "smoothed version" of the 
nth derivative of s(t) on various scales. 

In the presence of a singularity in the data at a particular point to, the scaling 
behavior of the wavelet coefficients is described by the Holder exponent a (to) 
as 

W s (a,t )~a a ^ , (5) 
in the limit a — > + . 

The continuous wavelet transform given in Eq. [2] is an extremely redundant 
representation, too costly for most practical applications. To characterize the 
singular behavior of functions, it is sufficient to consider the values and po- 
sition of the WTMM [14]. The latter is a point (a ,b ) on the scale-time 
plane, (a, b), where |W s (ao,&)| is locally maximum for b in the neighborhood 
of bo- These maxima are located along curves in the plane (a, b). However, 
the relationship (jSJ) in some cases is not appropriated to describe distribution 
functions with non-isolated singularities. The wavelet multifractal formalism 
may characterize fractal objects which cannot be completely described by a 
single fractal dimension. According to Bacry et al. [10J, an "optimal" partition 
function Z q (s,a) can be defined in terms of the WTMM. They consider the 
set of modulus maxima at a scale a as a covering of the singular support of 
s with wavelets of scale a. The partition function Z q measures the sum of all 
wavelet modulus maxima at a power q as follows 

Z q (s 1 a)=Y J \W s (a\(a))\\ (6) 

p 

where {b p (a)} p& z is the position of all local maxima of |W s (a, 6)| at a fixed 
scale a. This partition function is very close to the definition of the partition 
function described in [7j. It can be inferred from ([6]) that for q > the most 
pronounced modulus maxima will prevail, whereas for q < the lower ones 
will survive. For each g 6 1, the partition function is related to the scaling 
exponent r(q) in the following way Z q (s, a) ~ a T ^ 9 \ A linear behavior of r(q) 
indicates monofractality whereas nonlinear behavior indicates a multifractal 
signal. A fundamental result in the (wavelet) multifractal formalism states 
that the singularity (Holder) spectrum f(a) of the distribution function s(t) 
is the Legendre transform of r(q), i.e., 

a = r'(q), and f(a) = qa - r(q). (7) 

The Holder spectrum of dimensions f{ot) is a nonnegative convex function 
that is supported on the closed interval [a min , a max ], which is interpreted as the 
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Hausdorff fractal dimension of the subset of data characterized by the Holder 
exponent a [llj. The most "frequent" singularity, which corresponds to the 
maximum of f(a), occurs at the value of a(q = 0), whereas the boundary 
values of the support, a min for q > and a max for q < 0, correspond to the 
strongest and weakest singularity, respectively. 

The analyzing wavelets which are used most frequently are the successive 
derivatives of the Gaussian function 



because they are well localized both in space and frequency, and they remove 
the signal trends that can be approximated by polynomials up to the (n— l)th 
order [9]. In particular, our analysis was carried out with the Mexican hat 
wavelet ip^(t). 



4 Results and discussion 

We performed the multifractal wavelet analysis based on the WTMM method 
employing the biomass ANN estimates for two Agave spontaneous fermenta- 
tion processes. For the first one, we present the reconstructed biomass signal 
and the redox experimental signal in Fig. ([2]) (a,b), respectively. In general, 
for extracting the singularity structure of a signal according to wavelet trans- 
form techniques one should get the modulus maxima lines of the (continuous) 
wavelet transform. These are displayed in Fig. (j2J) (c,d) and one can imme- 
diately notice that for both reconstructed and real experimental redox data 
we do not obtain isolated singularities as those treated in our previous work 
PQ, on the contrary there is a high density of singularities. This fact induces 
us to think of the multifractal formalism which could be more appropriate in 
such cases. The complete evidence of the multifractal character of the signals 
is a scaling exponent r with two different slopes, where the point of bending 
in the r(q) spectrum is directly related to the most frequent singularity of 
the biomass data. It is worth noting that the locations of the maxima lines 
are almost the same for both the reconstructed biomass data and the direct 
experimental redox potential data, see the same figure (e). This indicates that 
the effect of the ANN procedure on the singularity structure of the redox 
signal is small and that we can be confident in the ANN-based total biomass 
estimates. Thus, assuming that these spontaneous phenomena are multifractal 
processes, we calculated in Fig. (j2J) (f,g) the two basic multifractal quantities, 
namely the r(q) and f{ct) spectra for both sets of data. We confirm the mul- 
tifractality of both signals since we indeed get a r spectrum with two slopes 
in both cases. The support of the singularity spectrum for the biomass signal 
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is a G [0.11, 1.36] with the maximum of the spectrum located at a pea k = 0.95 
that gives us the most frequent singularity fB(oc pea k) — 0.756, whereas in the 
case of the redox potential we get the support a G [0.1, 1.04] with a pea k = 0.984 
giving the most frequent singularity /i?(a pe afc) = 0.588. 

We repeated the same multifractal analysis for another data set of the Agave 
spontaneous fermentation process to examine the reproduction of our results 
of the previous case. The second analysis is displayed in Fig. (J3J). We again 
get the confirmation of the multifractal behaviour from the two-sloped r spec- 
trum, but of course as we deal with a different time series we have slightly 
different multifractal parameters. The support of the singularity spectrum 
for the biomass signal is now a G [0.188, 1.38] with the maximum of the 
spectrum located at a pea k = 0.84 providing as the most frequent singularity 
fB{.ctpeak) = 0.68, whereas in the case of the redox potential we get now the 
support a G [0.1,0.88] with a pea k = 0.7885 that gives for the most frequent 
singularity the value /j?(a pea fc) = 0.633. 

There is still one point to discuss, which is why the Holder support extends up 
to one for the experimental redox potential signal whereas it goes substantially 
beyond one in the case of the biomass signal in the two cases. We have found 
in the literature that a similar situation has been encountered by Humeau 
and collaborators [15] in the case of the laser Doppler flowmetry signals. The 
Holder exponents of their simulated data are greater than one as compared 
to their real signals for which they are less than one. They comment that in 
general the simulated data are smoother whereas the real signals are much 
more irregular. 



5 Conclusions 

Although mixtures of microorganisms may look natural in fermentation pro- 
cesses, the total biomass, which is the quantity most easy to deal with, does 
not provide definite evidence for the mixture by itself. However, in such pro- 
cesses the redox potential and the total biomass concentration are related 
fluctuating quantities thence the redox potential signal may contain relevant 
information on the microorganism metabolisms. This can be inferred through 
the multifractal techniques which has been specially designed for signals with 
significant irregular contributions. 

The main finding of this paper is that the spontaneous fermentation processes 
in the case of the Agave syrup can be characterized as a multifractal process 
and that the Holder exponents are a good measure of the degree of irregularity 
of the signals. Since the total biomass concentration and the redox potential 
signals are connected spontaneous time series we expect to get quite similar 
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Fig. 2. (a) Biomass signal for one of the six performed spontaneous fermentation 
processes of the Agave syrup (under the same experimental conditions) which was 
obtained by rebuilding from the redox potential data through the ANN procedure. 
The redox signal is displayed in (b). The maxima lines of the continuous wavelet 
transform of the biomass signal are shown in (c), of the transformed redox signal in 
(d), and in (e) one can see a comparison of their locations, (f) The corresponding 
r{q) spectrum, (g) The multifractal singularity spectrum /(a) of the biomass and 
redox signals of this Agave spontaneous fermentation process. 



singularity spectra for them, which indeed is the case. A detailed mapping of 
these Holder irregularities to the metabolic parameters awaits more experi- 
mental data processed through the multifractal procedure. In the literature 
one can find that many other natural processes have been approached by the 
multifractal means, such as the membrane ion activity of Kc a channels, i.e., 
their non-stationary dwell time series [16J, the atmospheric energy cascades 
[j~7] . and the point rainfall time series |18j . to name just a few. In the well- 
established multifractal techniques one can find embedded a useful amount of 
information and we strongly believe that the wavelet multifractal computa- 
tional scheme may provide an appropriate tool for quality control for this kind 
of spontaneous processes. The results of this paper together with our previous 
research on mixed growth [1] show that the wavelet methods are a very useful 
tool for detecting physical singularities in the fermentation processes. 
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Fig. 3. (a) Biomass signal for another of the six spontaneous fermentation processes 
of the Agave syrup that we performed by rebuilding it through the ANN from the 
redox potential data (same experimental conditions) presented in (b). The redox 
signal is displayed in (b) . The maxima lines of the continuous wavelet transform of 
the biomass signal are shown in (c) , of the transformed redox signal in (d) , and in (e) 
one can see a comparison of their locations, (f) The corresponding r(q) spectrum, 
(g) The multifractal singularity spectrum f{a) of the biomass and redox signals of 
this Agave spontaneous fermentation process. 
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